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ABSTRACT. We study a system of stochastically forced infinite-dimensional coupled harmonic oscillators. Although this 
system formally conserves energy and is not explicitly dissipative, we show that it has a nontrivial invariant probability mea- 
sure. This phenomenon, which has no finite dimensional equivalent, is due to the appearance of some anomalous dissipation 
mechanism which transports energy to infinity. This prevents the energy from building up locally and allows the system to 
converge to the invariant measure. The invariant measure is constructed explicitly and some of its properties are analyzed. 



1. Introduction 

Consider the infinite dimensional linear system of coupled stochastic differential equations 

f a n (t) = a n _!(i) - a n+ i(t) + <L iW^f) n e N 



<N ■ (l) 

>■ loo(t)=0 
. 

\Q \ where the forcing, W(t), is additive white noise and 6 n ,i denotes the Kronecker delta. Formally, the unforced system 

■ appears to conserve energy: 

o 

\Q ; (2) — I - a 2 n J = 2_j «nfln = a n (a n -i - a n+ i) = a n a n -x - a n a n+ i = 0. 

c"^ V nGN / nGN n£N nGN nGN 

Of course, this calculation is only formal as the rearrangement of the summations is justified only if the sequences 
<"*> f involved are absolutely convergent. Since the matrix which encodes the coupling in (HJ is real antisymmetric, any 

finite even dimensional truncation of this matrix has the same real Jordan canonical form as a system of uncoupled 
simple harmonic oscillators. 

Our goal in this paper is to address the following natural question: What is the long term behavior of (Q]l? More pre- 
cisely, does the forced system (Q~|i have an invariant measure or statistical steady state? Finite dimensional truncations 
might lead one to conjecture a negative answer to this question. This is due to the fact that, in the finite dimensional 
setting, there is no mechanism for energy dissipation. The infinite dimensional system, on the other hand, does indeed 
give rise to a nontrivial invariant probability measure which is not supported in 1%. Since there is no apparent dissipa- 
tive mechanism, we loosely refer to this phenomena as anomalous dissipation. The present authors have introduced a 
formal version of this notion in [3 1 and have used it to analyze a variety of models. 

The invariant probability measure described above is explicitly computable and is given by the formula: 



(3) a n {i) = \l-l G u n (t-s)dW(s), 

where 



G° n (s) = / U n ^(z)Vl-z^e 2lzs dz, 



U n {z) = i n U n {z), and U n {z) is the n th normalized Chebychev Polynomial of the second kind HHH. 

In this paper we show that (HJ does have a nontrivial invariant probability measure given by (0. We derive the 
explicit representation for this measure and compute its covariance structure. We also describe a natural class of initial 
data whose long term dynamics converge to (O. The analysis is technically elementary and involves only classical 
orthogonal polynomials and basic facts from stochastic calculus. The model analyzed here should be compared to the 
related models presented in |3 | whose structures are similar except that the coefficients of the a n _i and a„+i terms 
depend on n. Those models also exhibit anomalous dissipations, but the qualitative features of their invariant measures 

l 



differ from those of © (in particular, the models in [0 have a different covariance structure than 0). The analysis 
performed in |0, which is based on generating functions, is also different from the one for presented below. 



2. Invariant Measure 

We first analyze a damped version of (0 in the interest of guessing the inviscid invariant measure 0. Next, 
we prove that the damped invariant measure converges to the inviscid invariant measure 0. Consider the damped 
stochastically forced infinite dimensional system: 



a n (t) = a n _i(i) - a n+ i(t) ~ va n {t) + S nA W(t) 
_ oo (t) - 0, 

where v > is a parameter. We show that this equation has a stationary solution given by 

7T 

where 



(4) 



(5) 



G»(t-s)dW(s) 



G v n {s) 



U n -!{z)y/l- z 2 e-^- 2 " )s dz. 



Before establishing that (0 solves (0, we show that (0 is a well defined random variable. It is sufficient to prove 
that \G^(s)\ ds < oo. This is clear for the case v > 0. We present a short argument for the case v = 0. 
Let f%(z) = Un(z)y/l - z 2 e 2ilz where 7 <E [0,7r]. Note that \f^(z)\ 2 dz < 00 and independent of 7. Also 
note that G° (7 + kit) = J , f2_ 1 (z)e 2 '* lkz dz is the — k th Fourier coefficient of f2-i- The Plancharel theorem and 
compactness of [0, ir] imply that 

J2Kh + kn)\ 2 <C<™ 



for some constant C which is independent of 7. These observations give the desired estimate 



/ [Gi( s )] 2 ds <j2 K 
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(7 + kit) I c?7 



l/„ 7 -iW| 2 ^ 
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G°( 7 + fc7T)| 



c?7 



^ = ^ll/J-lll2 < Ctt < OO 



and show that (0 is a well defined random variable. In fact, equation (0 defines an infinite dimensional Gaussian 
Markov process, {<z^(t)}. The measure that induces is clearly invariant under time shifts of the driving Brownian 
motion. Thus, once we show that is a solution of equation (0 for all v > the previous observation implies that 
induces an invariant measure for (0 for all v > 0. Since the invariant Gaussian measure we constructed is explicit, we 
can calculate its covariance structure: This will be done in section [3] 

To show that 01 solves (0 for all v > 0, we use the following basic fact about Wiener integrals: For nice G, e.g. 
G S C^K) n #i(K+), g(t) = G(t - s) dW(s) satisfies 



(6) 



dg(t) = G(0) dW{t) 



G'{t-s)dW{s) 



dt. 



Denote by fi(dz) the measure \/l — z 2 dz. G^(0) = and the orthogonality of the {U n } in L 2 ([— 1, 1], /i) imply 
that G^(0) = for all n > 1. Calculating the G^'(t — s) and using the recurrence relation (0, we arrive at the relation 
-i 

{v - 2iz)U n -i(z)Vl - z 2 e- (v - 2lz){t - s) dz = G v n _ x {t - s) - G v n+1 {t - s) - vG v n {t - s). 



(7) G%{t-s) = - 

We recover the equations in (0 by applying (0 and to formula for the sequence {a^ } and using the recurrence 
relation 



(8) 



U n+ i(z) - U n -\{z) = 2izU„(z) . 
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Relation (0 follows directly from the well-known three term recurrence relation for the Chebychev polynomials: 

U n+1 (z) + U n -!(z) = 2zU n . 

This shows that (O does indeed provide a solution and stationary measure for the infinite dimensional coupled sys- 
tem ©. 

Finally, we explain the origin of formula <(3j for the {a^}. Notice that the coupling matrix of equation (Q]i is 
the Jacobi matrix associated to the three term recurrence of the Chebychev polynomials of the second kind J3]|4). 
Therefore, working with generating functions in the Chebychev polynomials will diagonalize ([T). With this in mind, 
consider 

oo 



Using the recurrence relation (0 and the facts Uo{z) = y § and U\{z) — 2izUo(z), one can calculate the time 
derivative of a v {z, t) as follows: 

oo 

Ct U (z, t) = ^ On(t)Un-l(z) 
n=l 
oc 

= E - - "«»(*) + KiW(tj\ Un-i(z) 

n=l 

oc oc 

= -va v {z, t) + J2 <-i(t)U n -i(z) - J2 <+i(t)Un-i{z) +U (z)W{t) 

n—1 n—1 

= -va v {z, t) + al(t)Ui(z) - 2iza v 1 (t)U (z) + 2iza v (z,t) +U (z)W(t) 
= -va v {z,t) + 2iza v {z,t) +U Q (z)W(t). 

This implies that a v (z, t) satisfies the stochastic differential equation 

(9) a v (z, t) = -{v + 2iz)a v (z, t) + J^W(t) 

with z viewed as a parameter varying in [—1,1]. The solution to the initial value problem for (O and t > s is 

(10) a v (z,t) = ip v (t- s,z)a"(z,s) + \j^- [ tp v (t - r,z)dW(r) 



where (p v (t, z) = exp ((2iz — v)t). Letting s — > — oo and assuming "nice initial conditions" and v > 0, we obtain 
the form of the invariant measure in formula ©. We note that although this is the way in which the expressions above 
were derived or "guessed", none of the results depend on this derivation. 



3. Covariance Structure of the Invariant Measure 

We compute the covariance structure of (fj), c°(m, n) — Ea^a°, by first computing the covariance structure of ©, 
c v (m, n) — Ea^a^, and taking the limit as v — > 0. The justification of this procedure requires two short steps: First, 
show that there is a sequence of Vk — > such that a^ fc (0) — > a° (0) almost surely as k — * oo; second, show that a v n is 
almost surely a uniformly continuous function of v E (0, M) for any M > 0. Therefore, the convergence of a y n — > a" 
as v — * holds almost surely. The following estimate and the Borel-Cantelli lemma complete the first step: 

(11) P (j<C(0) - a" (0)| > < k 2 \e-^ - l\ 2 \G° n (t)\ 2 dt. 

To apply the Borel-Cantelli lemma, simply choose so that the sum over k of the right hand side of (fTTT i is finite. 
For such a choice of Vk, a%° (0) — > a° (0) almost surely. To show that a^(0) is almost surely a uniformly continuous 
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function of v 6 (0, M), we appeal to the Kolmogorov continuity theorem and the following estimate. Fix < r) < |, 
p > z/ > 0, and observe that 



EK(0)-<(0)| 2 



e- 2l,t \\ - e- {p - v)t \ 2 \Gl(t)\ 2 dt 

e- 2vt \l-e~( p ~" )t \ 2 \G n {t)\ 2 dt+ / e- 2l,t \l-e- {p - u)t \ 2 \G° n {t)\ 2 dt 



l\p-v\-« 

< d|p - i/^-d + de-toW <C 3 \ P ~ v\ 2 ^\ 

where the last inequality holds for \p — v\ sufficiently small and G\,G%, and C% are constants which do not depend on 
rj, v, and p. This completes the justification for computing c°(m, n) by taking the limit of & '(to, n) as v — ► 0. 



E<(0)a-(0) 







(12) 



We now compute (n, n): 
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G^( S )G-(s)d S 



o 



W n _i(z)Vl-z 2 e-("- 2 ") a dz 



-l 



U n -i(z')Vl - z ,2 e-(»- 2lz ') s dz' 



-l 



f/.s 



l r i 



Z/ n _i (zK-i (z') Vl-zVl-z' 2 

1 A y/l - zVl - Z' 2 

U n -x(z)U n -i(z')— — —dzdz' 



e -2(v-c(z-z'))s ds 



dzdz' 



i-u-i 2v-2i{z-z') 
The fact that the dynamics of {a n } is real implies that one needs only to compute the real part of equation ( fT2l : 

.1 r i 



(13) 



C7„_i(z)C/„_ 1 (z') 



, Vl - zVl - z' 



U-i i/ 2 + (z-z') 2 



■dzdz'. 



sin(irO) 



Introducing the change of variables z = cos(7T0),z' = cos(7r#') and using the fact that U n (cos(Tr9)) = y^J sln ( 7r ("+ 1 ) e ) 
see for example [4|, leads to the following integral: 

1-1 smtnnz) smfirnz') sui(ttz) sintTTz') . 

dzdz . 



(14) d / (n,n)=2v 
Analyzing the limit v — ► 0, one finds that: 

(15) 



J Q V 2 + (C0S(7TZ) — C0S(7TZ')) 2 



c°(7i, n) = lim c v (n,n) — [ sin [ it— ) sin 2 [ 7r— ) 
*->o J \ 2/ V 2 / 



Note that linin—Kx, c°(rt, n) = — which immediately shows that the invariant measure is not supported on £2- We will 
discuss the implication of this fact in Section [4] For the moment, we only remark that the existence of this limit is 
consistent with the systems' invariance by translation except for the forcing and the boundary condition at n = 0. 
We compute the general covariance structure c v (m, n) in two steps: n — mis odd; n — m is even. First note that 



(16) E< (0)^(0) = 



(-1)* 



•n+m 



1 r l 



n ^{z)\j m ^{z')j\~z~ 2 ^J\ - {z'f \ +li ; z z 2 2 dzdz'. 

v l + (z — z'Y 



The cases when n — m is odd or even must be treated separately: 



Case 1: If n — m is odd, then 

(-!)« 



E<(0)<(0) 



E<C(0)<(0) 



(-1)' 



n+m+l 



1 r l 



I J -I 
1 rl 



U n ^(z)U m ^(z'Wl~z 2 ^l-{z'f 2 (Z . Z ' ] dzdz' 

v l + (z — z')^ 



E4-i(z)t/ m -i(z')Vl - z 2 ^l - (z') 2 



(*-*') 

v 2 + (z — z') 2 



dzdz' . 



Since n — m is odd, these two expressions have opposite signs. On the other hand, the fact that the dynamics of 
the {a n } is real implies that the two expressions must be equal. Therefore, c v (m 1 n) = which further implies that 

c°(m, n) = 0. 



Case 2: If n — m is even, then 
(17) E<(0)<(0) = 



(-l) n 

V >_p+m v 



U n -i(z)U m - 1 {z')y/l - zVl - i z ') 2 



v 2 + {z- z' f 



dzdz' . 



Once again, analyzing the limit of ( fT7] i as v 



c (to, n) 



lim c"(m, n) 

ls—>0 



leads to the formula for c°(m, n): 
r 2 

(-1 



(18) 



(-1)' 



1 



(n + rn) — 1 (ri — m) 2 



4. Basin of Attraction of the Invariant Measure 



Next, we prove that if the initial condition, {a„}, is in £ 2 , i.e. J^^Li l a «| 2 < 00 > tnen tne dynamics converges 
weakly to the invariant measure (f5]) for any i/ > 0: Any finite collection of coordinates, a v ix , a" , converges to ^ 
as the initial condition is pulled back to s = — oo. Assume that {a n } E £2 and construct the function 

00 

(19) a u Q {z) = ^a n Un-i{z). 

n=l 



If a v (z, s) = Q.q{z) is the initial condition at time s for equation ©, then the solution at time t > s is given by 
equation ( TTOb . Note that the dynamics is well defined for £2 initial data and the solution remains in £2 for all t < 00. 
To recover the solution at time t we simply use orthogonality of the {U n }- 

-1 



<(*) 



W„_i(z)a iy (z, t)yl — z 2 dz 



(20) 



U n -i{z)ip v (t — s, z)a v '(z, s)t/l — z 2 dz 



U n -i{z) 



ifi v (t-r, z)dW(r) 



yl - z 2 dz. 



Since U n ~\{z)a^(z)^/ \ — z 2 X[~i,i] ( z ) G standard Fourier analysis implies that the first integral vanishes as 

s — * —00. The second integral converges to the form of the invariant measure ©. The convergence is uniform if 
a finite collection of coordinates is fixed. Therefore, for any bounded cylinder function, convergence is established, 
which, in turn, establishes weak convergence of solutions with £2 initial data to the invariant measure Q. 

A drawback of the convergence result above is that, when v = 0, the invariant measure is not supported on £2- 
Hence, we may wonder about convergence of initial data which is not in £2 and is in the support of the invariant 
measure or even in 1^, This question, however, turns out to be quite complicated as the behavior of the initial value 
problem for (T} depends sensitively on the initial condition. Since (HJ is linear, it suffices to understand the solution of 
the unforced system with initial condition a n (s) = a°: 



(21) 



a n (t) = a„-i(i) 
a (t) = 0. 



a n +i(t), a„(0) 



?■/ G N 



The solution of the forced system with initial condition a n (s) = a° is then obtained by adding the solution of (HJ with 
dn{s) = to that of (I2TI 1. Cataloging the behavior of the solutions of (t2TT > is complicated, even if we restrict ourselves 
to initial condition in i^. The difficulty is immediately understandable if one notices that 

(22) a 2n = 1, a 2n +i =0, neN 

is a fixed point for (Bil l and belongs to i^. Similarly, one can find time-periodic solutions with arbitrary period T > 
which also belong to £00. For the sake of brevity, we will refrain from attempting a complete analysis of (fJTJ. We 
note, however, that such an analysis has been performed in detail in [3 1 for models related to (f2TT> . 
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